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We formulate the dynamical mean field theory directly in the continuum. For a given definition 
of the local Green's function, we show the existence of a unique functional, whose stationary point 
gives the physical local Green's function of the solid. We present the diagrammatic rules to calculate 
it perturbatively in the interaction. Inspired by the success of dynamical mean field calculations for 
model Hamiltonian systems, we present approximations to the exact saddle point equations which 
may be used in the strong coupling regime, by using mappings onto generalized quantum impurity 
models. 



I. INTRODUCTION 

The calculation of the electronic structure of solids starting from first principles has been the subject of intensive 
investigations in the past. The most widely used method is the density functional theory (DFT), used frequently in 
the local density approximation (for a recent review see Ref. |l|). In principle, this method is good for calculating 
ground state properties only. Time dependent extensions of density functional theoryo have also been formulated 
and allow the calculation of density and spin correlation response functions which contain information about excited 
state properties. This method has been very successful in understanding weakly correlated solids but has serious 
drawbacks when dealing with strongly correlated systems. A second general formulation|-pf the quantum many body 
problem is the Green's function method as formulated early on by Luttinger and Wardtl and Baym and KadanoflQ. 
In this method one sets up fully self consistent equations for the one particle Green's function using a functional i.e., 
an effective action for the Green's function, which has a well defined diagrammatic interpretation in terms of skeleton 
graphs. However, when considering a realistic band structure, this method becomes fairly intractable beyond the 
lowest order or GW approximationD, and is inadequate for strongly correlated problems. 

The purpose of this paper is to describe the foundation of a third approach, namely, the dynamical mean field 
method (DMFT) in the continuum, in relation to the two other techniques described above. Our main motivation 
stems from the success of DMFT in dealing with model Hamiltonians of strongly correlated electrons on a latticgj. In 
fact, there ase already several works attempting implementations of dynamical mean field ideas to electronic structure 
calculations^. This paper is an attempt to provide a theoretical underpinning to these computational efforts. 

The local spectral function is a central concept in dynamical mean field theory. The goal of the theory is to 
compute the one electron addition and removal spectra, and quantities which can be naturally expressed in terms of 
the local spectra. One can view the local spectral function as providing a frequency resolution of the local density. 
The hope is that by using a more complicated quantity to formulate the theory, one may be able to formulate fiexible 
approximations, which work well even in strongly correlated situations, where density functional theory in the local 
approximation is known to fail. Indeed, recent work on model calculations has shown that DMFT plays a very 
useful role in situations where significant changes in the local spectral function take place. These rearrangements 
of one particle spectral weight may not be accompanied by significant changes in the electron density, and may 
not be captured well in approximate implementations of density functional theory. This occurs in metal insulator 
transitions which are accompanied by very small symmetry changes or, small structural charges. Examples of this 
kind of transition are the paramagnetic insulator to paramagnetic metal transition in V2O3 and NiSei-xSx , the a 
to 7 transition in cerium, the trfuisitions form the /? to the 5 phase in plutonium or the more recently observed metal 
insulator transition in Yi-xTiQ^. We also expect the dynamical mean field theory to arovide useful descriptions of 
phases where local moments are well developed but not ordered such as in 5 plutoniumEl. 

In this paper, we will show that in general DMFT functionals are more cumbersome in their diagrammatic formu- 
lation than the Ward Luttinger functional for the full Greens function. Nevertheless, our motivation for presenting 
this formulation here is that DMFT is very promising in terms of practical implementation in very strongly correlated 
situations. The reason for that is that it can be facmulated in terms of impurity models embedded in a medium whose 
properties have to be evaluated self consistentlyli3. A large number of numerical methods for treating the resulting 
dynamical mean field equations have been developed in the context of model Hamiltonian calculationaj and these 
techniques can and should be extended to more realistic models of solids. 

Conceptually, the three techniques are very similar: one defines a functional F of the quantity of interest X (X being 
the density in the DFT case, the local Green's function (to be defined later) for DMFT or the full Green's function in 
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the Ward Luttinger technique) such that the minimization of T selects the physical value of the quantity of interest. 
In all the techniques, there is an exact functional Te^act having the required property but its explicit analytic form is 
very difficult to determine . In all the three cases, approximations to Texact are necessary to obtain practical results. 
We adopt a fairly general formulation to set up a diagrammatic expansion for the exact functional of the local Green's 
function, whose saddle point gives the physical local Green's function. These are presented in the first part of this 
paper and some possible approximate implementations are presented in the second part of this paper. 

Extending the more standard lattice model calculations, our work provides a foundation for doing ab-initio dynam- 
ical mean field theory plus electronic structure calculations in the continuum. The goal of our work is twofold. First 
we want to establish the exact diagrammatic content of F of the dynamical mean field formalism for realistic models 
of solids where electrons move in the periodic potential of a crystal, as a formal framework for performing electronic 
structure calculations or, in other words, to construct a functional for which a dynamical mean field approximation 
gives the exact answer to the electronic structure problem. This is the content of Sec. II. This perspective should 
be contrasted with an alternative formulationlil which will be presented in a separate publication, where dynamical 
mean field theory will be viewed as an approximate solution to the electronic structure calculation problem. Secondly, 
in Sec. Ill, encouraged by the success of dynamical mean field theory in the treatment of model hamiltoniansci, we 
suggest some useful approximations for this functional. This work should provide a conceptual framework for the 
recent efforts in combiningj-ealistic density functional calculations with techniques for solving impurity models in a 
self consistent environmentu. Finally, we notice that our work can be regarded as a generalization of Ref. |l^, where 
lattice DMFT was given a functional formulation in terms of a field which is conjugate to the local spectral function. 
This formulation has given useful insights into the physics of the Mott transition. 



II. EFFECTIVE ACTION FORMALISM 



In this section, we develop the effective action formalism for DMFT. This is in the same spirit as the effective action 
approach to Kohn-Sham density functional theory (where rtie bilinear operator is just— the electron density) which 
was developed in a series of publications by Fukuda it et alEJ and Valiev and Fernandotil. We first set the notation. 
The continuum is divided into unit cells labeled by R and r is the vector defined within a unit cell. A general "local" 
bilinear operator has the form 

0"^(r,r',T,r')= J dxdyK{r,r'\Ky)i;+{^,T)My^r') (1) 

where ip~^ and ip are the electron creation and annihilation operators and if is a kernel to be specified below, x, r are 
the space and imaginary time coordinates and a, (3 denote the spin indices. With the definition (TtO) = A{r, r', r, r') 
(|l|) leads to 

A"'3(r,r',r,r') = 1 dxrfy/f (r, r'|xy)G"^'(x, y, r, r') (2) 

In terms of the cell indices and unit vectors the space coordinates read x = R -I- u and y = R' + u'. Here, u, r', ru' 
are vectors inside a unit cell. The kernel K projects the full G onto a "local" Green's function A defined within a 
unit cell. By choosing the kernel K to be sufficiently localized around R — R' « 0, we extract or project out local 
information from the full Greens function. There is a great deal of flexibility in deciding what one would like to define 
as the "local" single particle Green's functions. The simplest choice would be 

K = Sn,n,5{r - r')6{r - u)5{r' - u') (3) 

which gives a resolution of the density with frequency 

'°(^) = ^ E E ^^P"""^ A""izu;„, r, r) (4) 

n a 

where A(iu;„, r, r) is the spectral function corresponding to A, uj„ are the Matsubara frequencies and /3 is the inverse 
temperature. Experience has suggested however, that it may be useful to retain local phase information among 
different orbitals in the local Greens functions. The definition 

K = <5R,R'(5(r - u)5(r' - u') (5) 
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gives us simply the restriction of the full Green's function to a given unit cell. For r = r', we regain the former 
choice of local spectral function. Note that (||) and hence the corresponding Green's function A is not invariant under 
simultaneous lattice translations of r and r' now viewed as coordinates defined throughout the whole crystal. A kernel 
with translation invariance is obtained by considering 

G,G' k 

Here the G are the reciprocal lattice vectors and k runs over the Brillouin zone. Rewriting the right hand side of (j^) 
in real space coordinates and summing over the reciprocal lattice vectors, we obtain a kernel invariant under lattice 
translations. 

X = Xexpik.(x-y)(5(r-u)(5(r' -u') (7) 

k 

where k runs over the Brillouin zone. Finally, given a set of tight binding orbitals, there is a natural kernel K, 
associated with this set. The full Green's function is expressed in a set of tight binding orbitals cfn as follows: 

G"'3(x,y)= J2 '/'i(x-Ri)</'m(y-R2)Gf^(Ri,R2) (8) 

hnR,iR,2 

Similarly, the local Green's function can be projected onto this set as 

A"P{r,r')^ X 0i(r-Ri),/.™(r'-Ri)Gf (Ri,Ri) (9) 

/mR,i 

Using can be rewritten as 

A"^(r,r')- J2 I '/'i(x-Ri)0i(r-Ri)./.™(y-R2)0,„(r'-R2)Gf^(x,y) (10) 

This leads to the following orbital decomposition for the kernel 

if(r,r'|xy)- ^ (/);(x - Ri)0;(r - Ri),/),„(y - R2)0™(r' - R2) (11) 

^mR,iR,2 

From these we see that different choices of the basis functions lead tcwlifferent kernels corresponding to different 
definitions of the "local" Green's function. We mention that recent worksQ incorporating DMFT concepts in electronic 
structure calculations have all used the LMTOla basis set. 

We now use the formalism developed inEj to construct the effective action F for the local Green's function defined 
by (^). These methods prescribe a scheme to construct an analog of the "Kohn-Sham" correlation potential calculated 
in the conventional density functional theoryEJ. The idea is to construct a perturbative expansion in the interaction 
strength, for the effective action F. We use the imaginary time formalism in the rest of the paper. The partition 
function Z for a system of interacting electrons coupled to a source field J via the operator O is 

Z[J"'3] =expM^[J"^] 

= j Di:D7p+eyi^[-S + j J"^(r, r', t, T')0"^(r, r', r, r')] (12) 

The coordinates x run over the entire range of the system and It is crucial to note that the indices af3 are not 
summed over in the second term in the exponential in (|T^). The quantity W[J] introduced above is the generator of 
all connected diagrams. S is the action of a system of electrons with Coulomb interactions between them. 

5 = j dxY,^t{x)[dr-^ + Vx{x)]i^c{x) (13) 
+ Y^J dxdx'^+ix)^+ix')Ucix - x')M^')Mx) (14) 

a;/3 
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In (|^, X = (x, r), Vx represents the ionic potential the electron moves in and the Coulomb interaction is Uc{x — x') 
|x — x'|^^(5(r — t'). The effective action F is defined by the Legendre transform 



T[A''f'] ^ W[Jap] - J"'^(r,/)A"'^(r,r') (15) 

r,r' 

A naive expectation might be that since A is just a "local" Green's function, the effective action functional F is given 
by the usual Baym-Kadanoff functional. Here we find that this is not the case, rather the restriction of r, r' to a 
unit cell induces non-trivial corrections to the F as will be illustrated below. To obtain F as a functional of A, the J 
dependent terms in ( |l5|) should be rewritten in terms of A. Since this is not tractable in the presence of interactions, 
we use the inversion methodtJ which involves an expansion of J, W in terms of the coupling constant of the interaction 
(for example, for coulomb interactions) and Jg, where Jq is the source in the absence of the interaction. The crucial 
ingredient of this method is that Jq is chosen so as to reproduce the correct local Green's function A"^ of the full 
interacting problem i.e., 



a0l 



a/3 



A"f^ (16) 



where Wq is the non- interacting free energy. Using (^6|), the series expansion for the and J[Jq^] can and 

rewritten in terms of A and thereby resulting in a series for the effective action F[yl"''] 

oo 

F[A"'3] = ^e2T,[A"^] (17) 

i=0 

The exact local function A of the interacting system is obtained by minimizing ( p^ ) 

-Ja^ - ^ (^^^ 



A. Diagrammatic rules for F 

The first step is the construction of the non- interacting " Kohn Sham" Green's function Gq [Jq] 



GfiR + r, R' + r') = 



dr-^ + vx{R + r)] W'5(r - r')<5a/3 - JS^{r, r')<5R,R.^ 



(19) 



Note that in the above expression R + r = (R + r, r^) and equivalently, r = r, r^. Jq^ is the Kohn-Sham hybridization 
function defined earlier. Jq is chosen such that Gq gives the prescribed local Greens function A of the full interacting 
problem. Using Gq, thej-Cjewalomb interaction lines e^Uc{R — R' + r — r') and a propagator D to be defined below and 
following the method ofHEJ, we can obtain F to all orders in the coupling e^. The zeroth order term is Wq = Tr log Gq 
and the first order term is given by the sum of the Hartree and Pock diagrams shown in Fig.l. The higher order 
terms which form the Kohn-Sham exchange-correlation potential in the conventional density functional theory are 
constructed according to the following diagrammatic rules. At every order : 

1. Draw all connected diagrams made of Kohn-Sham propagators Go and Coulomb interaction lines e^Uc{R—R'+r—r') 
with the corresponding weight factors. 

2. Eliminate all graphs which are one vertex irreducible i.e., separable by cutting a single Coulomb interaction line. 

3. For each two-particle reducible (2PR) graph (i.e. any graph that can be separated by cutting two propagator 
lines), perform the following operations.^ 

(a) Separate the graph by cutting the 2PR propagators. 

(b) For each of the two resulting graphs join two external propagators. 

(c) Connect the two graphs via the inverse two particle propagator D which is the inverse of the connected correlation 
function H 
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R,R',U,U' 

= -Y.^o''iR + r,U + u')Gt''iR + r',U + u) (20) 

R,U 

In the above expressions U + u = (U + u, t„). The inverse propagator D appears in the diagrammatics of the inversion 
method, as a resuh of corrections to Jq from the interaction which are constrained so that one reproduces the desired 
local Green's function A order by order in perturbation theory. 

(d) Repeat the procedure until no new graph is produced. 

(e) Sum up all the resulting graphs including the original graph. 




(a) (b) 

FIG. 1. a) Hartree and b) Fock diagrams 

In Fig. 2, we show the second order diagrams which contribute to F. Here the full lines denote particle propagators 
Go and the dashed line is the Coulomb vertex. Figs. 2a and 2b are examples of diagrams which are both one vertex 
and 2 particle irreducible. 




(b) (c) 

FIG. 2. Diagrams at second order 

However, since the diagram in Fig. 2c is two particle reducible, it transforms according to Rule 3 as shown in Fig. 3. 
The thick double line denotes the two particle inverse propagator D. 




R+r' R'+u' 

FIG. 3. Dressing of a two particle reducible diagram 



B. Self-consistency 



The self-consistent procedure for obtaining the full A is of the same kind as that prescribed for F[p] in Ref. |TJ. 

• Choose a Jo'^{r,r'). 

• Substitute this in ( p^ ) to obtain the Green's function Gq^{R + r,R' + r'). 

• Use this Gq^ to evaluate the series for F = Fq + + F^.^. Here Th denotes the Hartree contribution and Fj.^ 
represents the contribution of the Fock and all higher order diagrams to F. Determine the new "Kohn-Sham" 



Jq'^ defined as 



T"/3 _ 

Jg — 



6T 



H 



(21) 
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• Using the Jq determined by (^), go back to step 1 and iterate until self-consistency is achieved. 

In density functional theory, the second term in ( pT]) is called the exchange-correlation potential Vxc- The scheme pre- 
sented above, hinges on the fact that the Legendre transform defined in ( p^ is invertible. Mathematically, neglecting 
the spin indices, the transformation is invertible if 6^T/SA^ has no zero eigenvalues. Within the inversion method 
where F is evaluated perturbatively in the interaction strength, the invertibility condition described above reduces to 
the following equation : 

This yields (^6|) provided SJq/SA has no zero eigenvalues or equivalently, 



(5^2 



^ (23) 



Typically these transformations are invertible for generic values of the variable A. Non invertibility is related to 
a bifurcation of the stationary solution of the functional F and signals a phase transition. This was exploited in 
Ref. to analyze the energetics of a system near a Mott transition. The functional constructed here can be used 
in conjunction with the techniques of Ref. to provide more realistic descriptions of electronic phase transitions. 
Finally, we mention that we have not yet addressed rigorously, the questions corresponding to the issues of v and N 
representabilities in density functional theoryld i.e., what is the domain of definition of our functional or equivalently, 
what is the range of the function A7. The formal construction presented here is valid for A's which are Kohn Sham 
representable, and we know that the physical Green's function is Kohn Sham representable, with Jo being the physical 
"local" self energy. 



C. Example 

For the sake of clarity, we consider the case of spinless electrons with non-local interactions. In this case, 

n(r, r', u, u') ^ - Y,Go{R + r,U + u)Gq{R + /,[/ + u) (24) 

R,U 

We present the various terms upto second order in e^, contributing to T[A\ for the kernel of (^). The zeroth order 
term is Fo[A] = TrlogGo — Ji^A. At first order, we have the contributions from the Hartree and Fock diagrams shown 
in Figs, la and lb. 

Ff = /" drdr' Go{R + r,R + r)Go{R' + r', R' + r')Uc{R -R' + r-r') (25) 

R,R' 

The contribution from the Fock term is as follows: 



(26) 



Ff = /" drdr' Y Go{R + r, R' + r')Go{R' + r\R + r)Uc{R - R' + r - r') 

R,R' 

Similarly the diagrams which contribute to second order are shown in Figs.(2a-2c). 
V^^ j drdr'dudu' Y Go{R + r,Ri+u)Go{Ri+u,R + r) 

Uc{R -R' + r- r')Go{R' + r' , R'^ + u')Uc{Ri -R'^+u- u')Go{R[ + u' , R' + r') (27) 
J drdr'dudu' ^ Go(i? + r,Ri+ w)Go(i?i +u,R' + r')Uc{R -R' + r-r') 

Go(i?' + r', + u')Uc{Ri -R'l+u- u')Go{R'i +u',R + r) (28) 

(29) 
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^2 = J drdr'dudu' ^ Go{R + r,Ri + u)Go{Ri + u, R[ + u') 

Uc{R -R' + r- r')Go{R[ + u' , R! + r')Uc{Ri - R[ + u - u')Go{R' + r',R + r) (30) 

(31) 

Since Fig. 2c is 2PR, following the prescription (3a-3c) presented earlier, we subtract the following term from the 2PR 
diagram of ( |3l| ) 

^2 ^ j Go{Ri+ri,R\+i\)Go{Rx+rx,R + r)Go{R[+r\,R + r')Jl-\r,r' ,u,u') 

R,,R,' -R'l iR'"!^ if^2 iR'2 

Gq{R! + u,R2+ r2)Go(i?2 + r'2,R' + u')Go{R2 + r2, K + r^WdRi - R[ + n - +r[)Uc{R2 - i?2 + ra - +r^) (32) 

The integrations are over imaginary time and all the vectors within a unit cell. Therefore, the contribution to T at 
second order is 

r2-i(r^-rg) + (r^-rn (33) 

If we consider the case where the coordinates u, r scan the entire space i.e., K = 1 and A = G{x, y), we find that the 
structure of 11 yields F^'' = F2 , hence canceling the contribution of the 2PR diagram of Fig. 2c. Since such cancellations 
occur at all orders, the rules (1-3) imply that only two particle irreducible (2PI) are retained in the expansion for 



F. In addition, (16) implies that Go — G and Jq = G~^ — G~^, where G„ is the non-interacting Green's function in 
the absence of the source. Replacing Go by G in the above terms, we recover the Baym Kadanoff functional for the 
Green's function G 

F[G] = TrlogG - TrG^^G + Fi + /2[G] (34) 
where I2 is the sum of all connected two particle irreducible diagrams occurring at second and higher orders. 

D. generalization 

The formalism presented above can be used to construct the effective action functional for any arbitrary choice of 
the operator O or equivalently, the kernel K in (|l|). The rules (1-3) can be extended directly, provided K satisfies the 
following conditions: 

1. isr(r, r'|x,y) is invertible 

2. / X-i(r,r'|x,y)if(r,r'|w,z)drdr' = Jwx^y, 

3. / X"i(r,r'|x,y)if(u, u'|x,y)dxo?y = (Jruf^r'u' 

We can then use the rules (1-3) to construct T[A\ with the following replacements 

Go(i? + r,i?' + r') = [dr-^+vx{B. + T)^5nn'5{r-r')^ J dudu' J(u, u', r, T')i^(u, u'|R + r, R' + r') 

(35) 

and 

U{r,r',u,u') = J2 [ K{r,r'\Tl + a,Tl' + a.')K{u,u'\V + h,V' + h')(i:+{R + a)ij{R' + a')i:+iU + b)^{U' + b')) 

= Y [ K{r,r'\Tl + a,Tl' + a')K{u,u'\lJ + h,V' + h')Ga{R + a,U' + b')GoiR' + a',U + b) (36) 

The generalization to the spin dependent case is straightforward. 
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III. APPROXIMATIONS 



We now address the question of how to calculate T^c in As is the case in the usual density functional theory, 
the evaluation of all the diagrams that contribute to F and hence T^c is impossible and one has to take recourse to 
various approximation schemes depending on the strength of the interactions. 



A. Weak Coupling 



The simplest scheme would be to use truncated perturbation theory to second order, utilizing the expressions given 
in the previous section. Alternatively, one could use a random phase approximation. This approximation leads to a 
simplification of the rules for the effective action presented above. Since all RPA diagrams are by definition irreducible 
with respect to the coulomb vertex and two particle irreducible, we do not have to go through the complicated 
procedure outlined in rules (3a)- (3e) to deal with reducible diagrams. It would be interesting to investigate the 
equivalent of the local density approximation (LDA) for these new functionals, which corresponds to evaluating the 
diagrammatic expansion for T^c described in the previous section for a spectral function which is independent of the 
spatial coordinates i.e., spatially homogeneous. Unfortunately, the energy of the uniform electron gas expressed as a 
functional of its spectral function is not available and further work is required before a correspondence with LDA can 
be made. 



B. Strong Coupling - Impurity models 

The weak coupling approximations for F can be set up for an arbitrary kernel ii' in (|l|). On the other hand, different 
strong coupling approximations can be set up for differeHt choices of the kernel and are closely related to the formalism 
of dynamical mean field theory (DMFT) on the latticeO, where the local Greens function is represented as a Greens 
function of a generalized Anderson impurity model. In finite dimensions, the fact that the exact Green's function can 
be represented in terms of a local impurity model is an assumption. It suggests however, the implementation of the 
local Green's function defined by the kernel in (^). For simplicity, we consider the paramagnetic case in the following. 
Denoting the unit cell by B we represent the local Greens function defined by the kernel rtq) , as a generalized impurity 
model described by an action 



S{Go} = J drdr' j^dv j^dv'Y,^t{r)A-\r,r')i>^{r') 



IB J B 

+ fdrdr' [ dr [ dr' p{r)VsRir,r')p{r') (37) 

Jo J B JB 

Here r = (r,T) and the density p(r) — ?/'^(r)'i/;Q,(r). The interaction in the full continuum model is now replaced 
by an effective short ranged interaction Vsfl(r, r') = Vsr{YtY')6{t — r') within the cell B. For lattice models, such a 
mapping led one from the full Green's function Gij to a local Green's function Ga. In the continuum, the mapping 
leads to a similar "local" Green's function A{r, r') which is the Green's function defined within the unit cell. Referring 
back to Sec. II, we see that a simple generalization of DMFT ideas automatically leads us to consider ^ as the simplest 
choice for local Green's function in the continuum. Analogy with DMFT on a lattice, leads to write a self consistent 
equation for Ag by 

Ao\r, r') = [|: - 1^ - M - Vx{r) - (r)]5(T - r')<5(r - r') + A(r, /) (38) 

where Vh is the Hartree potential that arises naturally in such an impurity mapping and A(r, r') is the time-dependent 
"Weiss field of the problem". It describes the effect of all the other unit cells (which have been integrated out), on 
the selected unit cell. Note that this looks like a generalization of the multi-orbital impurity models if we assume that 
the orbital indices are the continuous variables r, r'. Using S{Ao) given in ( |37| ) we can compute the "local Green's 
function" A{r, r') of the unit cell defined by 

^'^'7isty^^'''°'--^(-'-^){^> (39) 



8 



Note that A depends on the Weiss field A. The self energy S[A] = — A~^{r,r') can be thought of as the sum 
of the Hartree term plus all the two particle irreducible skeleton graphs with interaction lines Vsij(r, r') and particle 
lines given by the "local" Green's function A{r,r'). Using the above, the self consistency condition for A(r, r') is 
obtained by requiring that 

A-\r, r')[A] = [-^ - ^ - - Vxir) - VH{r)[A]]S{r - r')^(r - r') + A(r, r') - I](r, r')[A] (40) 

This procedure of mapping the model onto an impurity model with interactions defined only within the unit cell 
is nothing but an approximation to the infinite series of diagrams in the pcrturbative expansion described earlier. 
Diagrammatically, the present method retains only local diagrams i.e., those whose internal propagators and vertices 
just connect points within the same unit cell. None of the diagrams include particle lines going from one unit cell to 
another. This is indeed reminiscent of the Hubbard model calculations in infinite dimensions where the self energy 
is completely local in space. A solution of this self-consistent equation for A is sufficient to specify the local A and 
hence other dynamical quantities. 

The |aJp|0¥e mapping can be generalized to a more general impurity model which goes under the name of extended 
DMFTIIjO. This captures some important aspects of the long range part of the Coulomb interactiontZl In the case 
of lattice models, one obtains a single impurity whose on-site interaction is dynamically and self-consistently screened 
due to the presence of long-ranged coulomb interactions in the original lattice model. This has been worked out 
explicitly for the lattice model in Ref. A straight forward generalization to the continuum is obtained by replacing 
the short-ranged interaction Vs/j(r, r') in (p7|), by the retarded time dependent interaction 

n-i(r, r') = VsR{r, r')6{T - r') + n{r, /) (41) 

In the same manner as above, the local density density correlation function 11 (p(r)p(r')) is given by 

n-\r,r')[n,A] = n-'-n[n,A] (42) 

where 11 is the irreducible part of 11. We, therefore, see that since 11 depends on the retarded correlations, the 
dynamical interaction Ho is also determined in a self-consistent manner. Moreover, in the presence of this retarded 
interaction, the self energy E = I][A, $7] i.e., it depends on A through Go and also on O through the interaction 
Hq. Therefore, ( ^0| ) and form a set of coupled equations for A and fl which have to be solved self-consistently. 
These are the generalizations of the extended DMFT equations obtained in Ref. [ij. The diagrammatic content of 
this approach is slightly different from that of the unscreened model since the vertex Vsr is now replaced by IIq. This 
implies that in the full continuum model, we retain all those "local" diagrams with internal particle lines A{r, r') and 
long ranged vertices which connect any pair of arbitrary points r -I- R and R' -f r'. 

To summarize, we see that DMFT ideas can be generalized to provide strong coupling approximations to the finite 
dimensional effective actions of the local Green's function A. We briefly mention that the above mappings can be 
generalized to include the effects of magnetic ordering. In this paper, we have presented the generalized impurity 
model only for the simplest kernel and it would be interesting to study whether there exist other generalized impurity 
models for arbitrary choices of the local Green's functions. In particular we notice that current DMFT calculations 
using realistic multiband HamiltoniansQ based on LMTO orbitals can be viewed as approximate realizations of the 
functional of the local Green's function based on the choice of the kernel (^l]), where the tight binding basis defining 
the kernel are currently obtained from the self consistent solutions of a density functional calculation. 

IV. CONCLUSIONS 

In this paper, we have described a means of studying electronic structure problems defined in the continuum, 
involving effective action formalisms and the principles of dynamical mean field theories. We first outlined the 
method to formulate the effective action functional for a general "local" Green's function. These local functions are 
good probes of the dynamics of the system including metal insulator transitions. Though we present a scheme to 
calculate these functionals to every order in the interaction, as is the case with standard density functional theory, 
an exact evaluation of these functionals is practically intractable. We have, therefore, described various weak and 
strong coupling approximations. The strong coupling approximations are a generalization of the DMFT ideas on 
a lattice to the continuum. Recently, some of these-generalized DMFT ideas were applied to the problem of the 
Hubbard model with coulomb long range interactions O. We found that the long range interactions screen the effective 
onsite interaction in such a manner as to make the continuous Mott transition seen in that system, discontinuous. 
These results suggest that a direct implementation of DMFT to an electronic structure problem can produce results 
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which are not just quantitatively but also qualitatively different from those obtained from conventional DMFT model 
calculations. 
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